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As an application of quantum fluid mechanics, we consider the drag force exerted on a sphere by 
an ultra-dilute gas. Quantum mechanical diffraction scattering theory enters in that regime wherein 
the mean free path of a molecule in the gas is large compared with the sphere radius. The drag force 
is computed in a model specified by the "sticking fraction" of events in which a gaseous molecule 
is adsorbed by the spherical surface. Classical inelastic scattering theory is shown to be inadequate 
for physically reasonable sticking fraction values. The quantum mechanical scattering drag force is 
exhibited theoretically and compared with experimental data. 
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I. INTRODUCTION 

Quantum fluid mechanical effects are very often con- 
sidered to be negligible. Except for some very special 
cases 0, such as the study of the superfluid phases of 
helium]^, , quantum fluid mechanics is rarely consid- 
ered. Here we wish to consider an important exception 
to this rule; i.e. the drag force exerted on a moving sphere 
by a highly rarified gas. We wish to consider the case in 
which the mean free path of a gas molecule is large on the 
length scale of the sphere radius. For example, a very rar- 
efied gas0 exists above the upper atmosphere. Meteors 
or spaceships on first entering such an atmosphere 0, 
may approximate the situation to be studied in this work. 

A rarified gas will exert a drag force on a moving 
sphere. If the mean free path of a molecule in the gas is 
small compared with the radius of the sphere, then the 
drag is due to the viscosity of the gas. If the gas is further 
diluted so that the mean free path of a molecule is much 
larger than the sphere radius, then the drag force in a 
kinetic theory picture depends on the notion of a stick- 
ing fraction /, i.e. the fraction < / < f of molecules 
incident on a surface which stick and thermalize to the 
temperature of the sphere before later evaporating. We 
shall later presume that those molecules which do not 
stick to the surface are specularly reflected. The central 
result of such a Knudsen model 6] is a relationship be- 
tween the sticking fraction and the slip drag force on the 
sphere. 

In past treatments of the scattering of molecules off 
the sphere, the classical scattering theory has been 
employed [9j. A central result of our work is that the 
classical scattering theory is inadequate and should be 
replaced by quantum scattering theory. The drag coef- 
ficient is proportional to the transport cross section for 
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molecules to scatter off the sphere. Due to diffraction ef- 
fects, the quantum mechanical cross section will differ ap- 
preciably from the classical cross section. The incoming 
molecules are sufficiently large and fast for their quantum 
wavelengths to be very small on the scale of the sphere 
radius. This is a necessary, but not sufficient condition 
for the classical limit to be obtained. The classical limit 
still fails to hold true due to diffractive shadow scattering 
which persists to even the smallest wavelengths. Diffrac- 
tion effects increase the total cross section by roughly a 
factor two. The factor is exactly two for purely specu- 
lar reflection. Millikan^(j made a series of experiments 
measuring the drag force on oil droplets by rarefied gases 
and found an effective transport cross section given by 
a m ~ 1.377TO wherein a is the sphere radius. It is not 
very easy to understand why such an experimental cross 
section is larger than the standard classical geometric 
cross section ira 2 . Unfortunately we do not know of any 
recent experiments, which attempt to measure the drag 
forces on spheres in a similar Knudsen regime ^J> We 
thereby use the reliable Millikan experimental results to 
compare with theory. 

In Seem firstly we review the classical kinetic the- 
ory of the drag coefficient on a sphere for the case of 
completely elastic classical specular reflection. Secondly, 
the completely absorptive limit is briefly discussed. In 
Sec lIIII we firstly review the kinetic theory of the drag 
coefficient on a sphere for the case of completely elastic 
quantum specular reflection. The quantum kinetic theory 
for the drag coefficient will be calculated. Secondly, the 
completely quantum absorptive limit is briefly discussed. 
In Sec lIVI a general sticking fraction model will be dis- 
cussed. The results obtained thereof will be compared to 
experiment. In the concluding Sec0 the implication of 
these results for the picture of quantum turbulent back- 
flow will be discussed. Quantum diffraction effects lead 
to a strong forward peak in the differential cross section. 
In quantum hydrodynamic terms, the forward scattering 
peak translates into a thin quantum trail of fluid which 
will flow behind a moving sphere. The trail is to those 
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particles which undergo specular reflection from the sur- 
face whose probability exhibits strong diffraction and in- 
terference effects. 



II. CLASSICAL DRAG FORCE ON A SPHERE 

Many of the kinetic models applied to drag forces in 
ultra-rarefied gases involve the notion of a sticking frac- 
tion / defined as follows: (i) (1 — /) is the probability 
that a molecule will elastically scatter off the surface in a 
specular fashion, (ii) / is the probability that a molecule 
will stick to the surface and thermalize to the substrate 
temperature before finally being evaporated back into the 
gas. Below, the drag force on a sphere of radius a will 
be reviewed using purely classical scattering theory and 
classical kinetic theory. Two opposite limits will be dis- 
cussed; namely (A) the purely elastic specular reflection 
limit, i.e. / — 0, and (B) the purely inelastic absorption 
limit, i.e. / = 1. 



A. Classical Elastic Specular Reflections 

Consider a particle moving in the gas in the direction 
of a unit vector u. If the particle hits the sphere and spec- 
ularly scatters through an angle 9, then the component 
of the momentum transfer along the original direction is 
given by 



u- (Pi-P/) = |p|(l -cos0). 



(1) 



the transport cross section is the same as total cross 
section([l4|); i.e. 



<7i = a = na 



(classical elastic hard sphere). (6) 



The drag force for the above cross sections is [l5L IIgIITtI 
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where 
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and erf(s) is the error function |2lj. When the velocity of 
the object is low, then the drag force obeysj^ 
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(10) 



The mean speed of a gas particle is denoted by c. 



If one now has n gas particles of mass m per unit vol- 
ume with a distribution of momenta p = pu, then the 
momentum transfers give rise to a drag force 



F = ( np^- / (1 - cos9)da 



(2) 



In Eq.(J2J), da is the differential cross section of scattering 
from the sphere into a solid angle dVt about the angle 0. 
The brackets (. . .) denote averaging over all momenta p 
with a Maxwellian distribution corresponding to a mean 
gas "wind velocity" v. The transport cross section is 
usually defined 13] as 



(1 — cos9)da, 



so that drag force may be written 



P 
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(3) 
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For classical specular reflection from a spherical sur- 
face, 



da 
dtt 



specular 



= — (classical hard sphere), (5) 



B. Classical Purely Inelastic Absorption 

By purely inelastic classical absorption we mean that 
any incoming particle whose impact parameter is less 
than the sphere radius sticks to the spherical surface with 
probability one; i.e. f — 1. The particle may much later 
be re-emitted after thermal equilibrium with the sphere 
is established. This kind of evaporation implies the re- 
emission of gas particles with a Maxwellian distribution 
with zero mean velocity in the reference frame of the 
sphere. The spherical nature of the re-emission implies 
the equality of the transport cross section in Eq.Q and 
total cross section. Specifically, 



ira classical. 



(11) 



It then follows that the drag force on a sphere due to 
purely classical elastic collisions coincides with the drag 
force due to purely classical inelastic collisions 22] . The 
drag force is again, 



P 1 2/~t 2 

F = -pna C F v , 
wherein Cf is defined in Eq.JSJl. 



(12) 
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III. QUANTUM DRAG FORCE ON A SPHERE 

The elastic amplitude for a gas molecule to scatter off 
a sphere may be expanded into partial waves 



Mp) 



2ip 



£(2l + l)[S,(p)-l]i5(cos0). (13) 



The total cross section can be decomposed into an elastic 
plus an inelastic part 

cr to t(p) = <Jei{p) + <y in (p). (14) 
The total cross section follows from the "optical theorem" 
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The elastic cross section is determined by 

da e . 



= \Mp)Y 
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^(2? + l)|l-5 i (p)| 2 , (16) 
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so that 



^n(p) =[^t) £(2* + !)[!- l 5 l(P)l 3 ]' 

^ P '1=0 

Thus, the probability wf of elastic scattering and the 
probability wf 1 of inelastic scattering in a given partial 
wave are determined, respectively, by 

wf = \Si(p)\ 2 and wi n = l-\S l (p)\ 2 . (18) 



A. Quantum Pure Elastic Scattering 

If the probability of elastic scattering in a partial wave 
is unity, then one may define phase shifts {5i(p)} via 



•jf = 1 implies Slip) = e 2l5l{p) . 



(19) 



The elastic cross section Eo. l|16|) yields the quantum re- 
sult 
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^(2Z + l)sin 2 5,(;p). (20) 



z=o 



The quantum mechanical effect of the drag force is de- 
termined by the transport scattering cross section <Ji(p). 
Eqs.Q, JXtj! and JTSJ implvpt 
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^/sin 2 ^^)-^)]. (21) 



The hard sphere phase shifts[24| are 

jiipa/h) 



tan 8i (p) 



ni(pa/h) ' 



(22) 



jl and ni are, respectively, the spherical Bessel and Neu- 
mann functions. 

In the high energy limit pa » h: Firstly, the asymp- 
totic form of the phase shift is given by 



5i{p) — > {—pa/Ti) + (Itt/2) as p — > oo. 



(23) 



Secondly, the partial wave summation cuts off at Til sa pa, 
where a is the radius of the sphere and p is the particle 
momentum. By inserting Eq. (|23|l into Eq. l|21|l one ob- 
tains 



lim <7i(p) = 2na 2 



(24) 



Comparing the classical specular reflection transport co- 
efficient in Eq.© with the quantum specular reflection 
Ea.Q24|) one finds the drag force ratio 



(Ji (quantum specular) 
cti (classical specular) 
_F(quantum specular) 
F (classical specular) 



(25) 



The argument for the famous factor of two [23 between 
the quantum and classical cross sections is that there 
exists an interference of amplitudes between the scattered 
wave and the incoming wave. This interference creates 
a peak in the forward direction. This effect is closely 
analogous to Fresnel diffraction in optics, wherein the 
limit to geometric optics cannot really be achieved. Our 
point here is that the so called classical limit h — > 
cannot really be achieved because of diffraction effects 
when the particle scatters off the sphere. From Eq.QJ 
one gets the final Eq. i|25[l . The drag force on a sphere 
due to elastic scattering is twice as large as the classical 
value if quantum mechanical diffraction effects are taken 
into account. 



B. Quantum Pure Absorptive Scattering 

Pure absorptive scattering takes place when the inelas- 
tic cross section is at maximum. Ens.i|17jl and (|18f) imply 
that pure absorptive scattering in the I th partial wave 
occurs^ when Si = 0. If all the partial waves scat- 
ter in a purely absorptive manner, then Sy^a) ~ and 



S, 



(l>ka) 



1. The total cross section becomes a to t ~ 2na 



The elastic cross section is thereby equal to the inelastic 
cross section; i.e. 



and atot — 2ixa 



(26) 



1=0 



As in the classical result, the inelastic cross section is 
equal to the inelastic transport cross section. The ratio 
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between the force due to quantum scattering and classical 
scattering is 



o\ (quantum absorptive) 
a\ (classical absorptive) 
F (quantum absorptive) 
^(classical absorptive) 



= 2, 



(27) 



The drag force due to quantum purely inelastic scattering 
is twice as large as the classical drag force. 



IV. STICKING FRACTION MODEL 

For purely elastic scattering, the 5-matrix eigenvalues 
obey the unitary condition \SA 2 = 1. In terms of the 
Heitler AT-matrix eigenvalues [221 13, 



Si = 



1 — iTrKi 
1 + IttKi 



(28) 



the unitary condition is enforced by requiring that K\ be 
real. In the most simple sticking fraction model, the in- 
clastic processes are described with imaginary if-matrix 
eigenvalues, 



txKi = —irji where rji > 0. 



(29) 



In the I partial wave, a fraction wf 1 of incident particles 
"stick" to the sphere and a fraction wf specularly reflect 
from the sphere. In detail, Eqs.JIBJ), GHl and (J2HJ) imply 



>"i 



(i-m) 2 

4t# 



(non — sticking) , 
(sticking). (30) 



Finally, we presume a single value for the sticking fraction 

wf 1 « / if hi < pa, 

w\ n w if U> pa. (31) 

Eas. ljTBjl - ijTSjl and P ^ -ipT )! imply the central results of 
the simple sticking fraction model; i.e. 



(ira 2 
(™ 2 )/ 



) i-Vw 



(32) 



The physical significance of the central quantum scatter- 
ing Ea. H32|l is as follows: (i) The cross section for the inci- 
dent particle to stick to the sphere is simply the classical 
cross section times sticking probability, i.e. a m = (to 2 )/ 
which would also be valid in a classical scattering context, 
(ii) The elastic cross section a el — (na 2 ) [l — \/l — /] 
is in part due to diffractive scattering which is a purely 
quantum mechanical process having no classical counter- 
part, (iii) The total cross section in the sticking fraction 
model is 



2-na 1 



1 — y/l — / (Quantum Sticking). (33) 



On the other hand, for a classical sticking model, the 
total cross section is geometrical; 



classical 



na 2 (Classical Sticking). 



(34) 



The physical kinetics of sticking and specular reflection 
are such that the mean momentum transferred to the 
sphere are equal for elastic and inelastic events The fric- 
tion drag coefficient is thereby 



lim 



r 



'1W \ atot 
9nm / 



(35) 



The quantum expression for the drag coefficient in the 
sticking model is thereby 



(27TO 2 



128k B T 
97rm 



(36) 



The sticking fraction for very dilute air molecules 
bouncing off an oil surface had been measured by two 
different methods: (i) In the rolling cylinder methodp9|. 
the drag force on one rotating cylinder due to the dilute 
gas between it and another nearby concentric stationary 
cylinder is employed to measure the sticking fraction /. 
(ii) In the falling drop method [To|. the drag coefficient 
r on a falling spherical oil droplet is measured. T then 
determines the sticking fraction /. The droplet should 
have a small radius compared with the mean free path 
length of an atom in the gas. Agreement between these 
two methods is obtained if the fully quantum mechani- 
cal scattering theoretical Ea. (|36|) is used in the analysis 
of the data. The experimental results are shown in the 
table below: 



Experimental Method 


Sticking Fraction / 


Rotating Cylinder 
Falling Droplet 


0.893 
0.901 



The close agreement corresponds to a measured total 
cross section for a droplet in which the Knudsen num- 
ber Kn ^> 1. Recall the conventional definition!^ of this 
number, 



Kn=± 
a 



(37) 



where A is the mean free path of a molecule in the rarified 
gas. The experimental number deduced from Ea. (|35|) is 



1.377T0f 



(38) 



from which / may be deduced from Eci. (|33|l . The experi- 
mental cross section a tot is larger than the classical total 
cross section cr'° ( ^ ;ssica j = ira 2 due to quantum diffrac- 
tion effects. The agreement between the experimental 
methods is not entirely satisfactory since only for one oil 
droplet is data available to us in which Kn w 10 2 ^> 1. 
However, the data that does exist seems to demand quan- 
tum diffraction effects for the scattering of a molecule off 
the sphere. Further experiments would be of interest for 
probing the accuracy of these results. 
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V. CONCLUSION 

The drag force on a sphere moving through a highly 
rarefied gas with large Knudsen number has been dis- 
cussed. It was found that quantum mechanics substan- 
tially alters the drag force. Quantum mechanics enters 
into the computation via diffraction effects in the cross 
section for molecules scattering from the sphere. In the 
extreme case, it was found that if the molecules of the 
gas are completely elastic or completely inelastic, then 
the cross section is twice that found if only classical me- 
chanics is taken into account. On the other hand, when 
there are both elastic and inelastic processes then the 
total cross section may be only somewhat higher than 
the classical result. When collisions are mostly inelastic 



and absorptive, the cross sections are mainly determined 
by the sticking fraction / for the molecule to thcrmalizc 
on the oil drop surface. Whereas classical elastic scat- 
tering is roughly isotropic, leaving an "empty shadow" 
behind the sphere, the elastic quantum cross section ex- 
hibits diffraction patterns strongly peaking in the forward 
direction of the gas flow. A quantum wake with a narrow 
stream of particles will thus appear behind the sphere. 

By independent measurements of the sticking fraction 
and the drag force, the quantum mechanical theory has 
been experimentally shown to be a more accurate de- 
scription of the flow than the classical cross section. How- 
ever, experimental data available to us in the Kn 3> 1 
regime is somewhat limited. Further experiments would 
be of great importance. 
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